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We calculate the in-plane Seebeck coefhcient of Sr 2 Ru 04 within a framework combining electronic 
structure and dynamical mean-field theory. We show that its temperature-dependence is consistent 
with entropic considerations embodied in the Kelvin formula, and that it provides a meaningful 
probe of the crossover out of the Fermi liquid regime into an incoherent metal. This crossover 
proceeds in two stages: the entropy of spin degrees of freedom is released around room-temperature 
while orbital degrees of freedom remain quenched up to much higher temperatures. This is confirmed 
by a direct calculation of the corresponding susceptibilities, and is a hallmark of ‘Hund’s metals’. 
We also calculate the c-axis thermopower, and predict that it exceeds substantially the in-plane 
one at high-temperature, a peculiar behaviour which originates from an interlayer ’hole-filtering’ 
mechanism. 


When a thermal gradient VT is established in a ma¬ 
terial, an electric field —VV is also generated. This 
thermoelectric effect can be used for solid-state refrig¬ 
eration and waste-heat recovery. The Seebeck coefficient 
(or thermopower) is the ratio a = —W/VT, measured 
under the condition that no electrical current flows. It 
not only determines the suitability of the material for 
thermoelectric applications, but is also a useful funda¬ 
mental characterization of its electronic state. In metals 
with strong electronic correlations, the slope of the ther¬ 
mopower at low temperature has been shown to scale 
with the linear coefficient of the specific heat [1]. This is 
an example of often noted but poorly understood relation 
between the thermopower and the entropy S. A precise 
relation has been established only for a free electron gas 
at low temperature, in which case, a = —Sjne [I] , and 
in the high-temperature atomic limit where the Heikes 
formula au = (fn)_E [1: 2]. This formula 

has been used to interpret the saturation of the ther¬ 
mopower at high-temperature in several transition-metal 
oxides [3, 5, 6]. More recently, it was suggested [7, 8] that 
the ‘Kelvin formula’ ax = (ff)y is a good approxi¬ 
mation to the thermopower. 

Here we consider ruthenates, especially Sr2Ru04 . We 
show that the thermopower provides key insights into 
the degrees of freedom which are relevant to the physics 
of these materials, and clarify when and how entropic 
considerations apply. Sr2Ru04 is arguably the cleanest 
and best documented among strongly correlated oxides, 
with extensive studies of its low-temperature Fermi liq¬ 
uid (FL) behavior [9] and unconventional superconduc¬ 
tivity [10]. Recent work shifted the focus to proper¬ 
ties at higher temperature and energy. It was shown 
that the low value of the scale Tfl ^ 25 K [11] be¬ 
low which Fermi liquid behaviour is observed, and the 
large effective mass enhancement, are due to the Hund’s 
rule coupling [5]. Ruthenates belong to a broader class 
of compounds that notably includes iron pnictides and 


have been called ‘Hund’s metals’ [13-15] (for a review, 
see [16]). Quasiparticle excitations, as revealed by a well- 
defined spectral function peak, were shown to persist well 
above Tpu , and the changes in the dispersion of these re¬ 
silient quasiparticles [17] as a function of temperature [18] 
and frequency [19] explain the deviations from Fermi liq¬ 
uid behavior. Overall, these results suggest that one can 
characterize the crossover into the incoherent regime bet¬ 
ter than previously imaginable. 

The ab-plane thermopower of ruthenates initialy in¬ 
creases linearly with T, as expected in a FL, and sat¬ 
urates for T ~ 300 K at a value of order 25-35 ^V/K 
(see Fig. 1 for Sr2Ru04and SrRuOa). As noted in 
Refs. [5, 20, 21], this value depends weakly on the cation, 
the lattice structure and the doping. This universal¬ 
ity is intriguing: does the value of thermopower reflect 
the active degrees of freedom in the metallic state above 
Tfl and thereby reveal its physical nature ? Boldly ne¬ 
glecting configurational and orbital contributions to the 
Heikes formula, the authors of Ref. [5] proposed that only 
spin degrees of freedom must be retained in the atomic 
entropy. May atomic considerations apply to an itinerant 
metal, and if so what happened to the neglected degrees 
of freedom ? 

In this letter, we use electronic structure and dy¬ 
namical mean-field theory (LDA-I-DMFT) to perform 
materials-realistic calculations of the thermopower of 
Sr2Ru04 based on the Kubo transport formalism. We 
show that the observed value of the room-temperature 
thermopower is explained by the fact that the spin de¬ 
grees of freedom are fluctuating in this range of temper¬ 
ature, while the orbital moments remain quenched up to 
~ 500 K. At this temperature, the thermopower displays 
a downturn, indicating that the decoherence proceeds via 
two distinct stages, a hallmark of Hund’s metals. We find 
that the in-plane thermopower is well approximated by 
the Kelvin formula for all temperatures, establishing a 
connection to entropic considerations. We also calculate 
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Figure 1. Temperature dependence of the in-plane ther¬ 
mopower in Sr 2 Ru 04 . LDA-I-DMFT Kubo results are com¬ 
pared to experiment at low T [22, 23] (full circles) and at 
high-T [24] (open circles), and to the result of the LDA- 
Boltzmann theory (crosses). The entropic Kelvin estimate 
uk = —l/edS/dn\T (plain line) is close to the full Kubo re¬ 
sult. The dashed horizontal line indicates the Heikes-like es¬ 
timate derived in the text including only spin degrees of free¬ 
dom. A similar behaviour is found experimentally for other 
ruthenates (triangles: data from Ref. [20] for SrRuOa). 


the c-axis thermopower and predict a strong enhance¬ 
ment at high temperatures, which we explain by a hole- 
filtering mechanism resulting from the crystal structure 
of Sr 2 Ru 04 . 

We use the LDA-I-DMFT method, in the implemen¬ 
tation of Refs. [4, 25, 27], as applied to Sr 2 Ru 04 in 
earlier work [5, 19]. The rotationally invariant inter¬ 
action H = {U - 3J)NiN - l)/2 - 2JS^ - Jp 12 
was applied to the t 2 g atomic shell, with S and 
L the total spin and orbital pseudo-spin opera¬ 
tors, respectively. The same interaction parameters 
U = 2.3eV and J = 0.4eV as in Refs. [5, 19] 
were used. The thermopower in DMFT is a = 
—fcs/e / dujT{uj)/3uj{—df /duj)/ f dujT{uj){—df/duj), 
where f{uj) is the Fermi function, /3 = l/ksT, and the 
transport function reads: 


Tp) = ^ ^ Tr vj:A^{uj)vj:A^{uj) 


V 


( 1 ) 


Here, H is a normalization volume, A^(oj) the spectral 
function matrix, and vj: the band velocities in the x- 
direction (in-plane response) or z-direction (out of plane). 
Self-energies obtained by a continuous-time Monte-Carlo 
impurity solver [27, 28] were analytically continued us¬ 
ing both a stochastic maximum entropy method [29] and 
Fade approximants [30]. The total energies were calcu¬ 
lated from charge self-consistent calculations as imple¬ 
mented in Ref. [4]. 

The calculated thermopower is shown on Fig. 1 (full 
squares with errorbars). It increases linearly with tem¬ 


perature, reaches a maximum close to 500 K and then 
slowly diminishes as the temperature is further increased. 
Overall, the theoretical values agree reasonably well with 
experimental data.Above room-temperature, our results 
exceed the experimental values, even when the systematic 
error due to analytical continuation is taken into account. 
(The estimation of the errorbars that are indicated on 
Fig. 1 is discussed in the supplemental material [30] in 
which the overall consistency with Fade continuation is 
also presented). Also displayed on Fig. 1 is the LDA- 
Boltzmann estimate [31]. As expected, the slope of a{T) 
at low-T is then too small by a factor of about four, con¬ 
sistent with the mass renormalization due to correlations. 

We now consider the Kelvin approximation uk = 
that arises [7] if the slow (thermo¬ 
dynamic) instead of the fast (transport) limit is used 
in the evaluation of Onsager coefficients. Alternatively, 
it can be derived using a physical argument, requesting 
that the density gradient vanishes instead of the current 
i.e. Vn(/x,T) = |^|tV/ 4-|- |y]^Vr = 0, and using the 
Maxwell relation “ V/x/eVT. 

Note that in the high-T limit where /i(T) oc T, the 
Kelvin expression coincides with the Heikes formula 
an = k'/dT, but that it has a greater degree of gener¬ 
ality and can be investigated at any temperature. As 
shown on Fig 1, the Kelvin expression (obtained from 
a numerical derivative of the LDA-I-DMFT chemical po¬ 
tential [32]) agrees with the Kubo result remarkably well. 
This finding extends to a real material the good agree¬ 
ment previously noted in model calculations [7, 17, 33]. 

The success of the Kelvin approximation hints at an 
entropic interpretation of the thermopower. In order to 
identify which degrees of freedom are active at a given 
temperature, we calculated the local spin and orbital sus¬ 
ceptibilities, displayed on Fig. 2 (a,b) as a function of 
temperature and frequency. The spin susceptibility has 
a Curie-like behaviour, indicating fluctuating spin mo¬ 
ments which are quenched only below the very low scale 
Ffl [5], which can be interpreted in DMFT as the Kondo 
temperature associated with the atomic t 2 g shell coupled 
to its environment. Orbital susceptibilities behave in 
a drastically different manner. They are much smaller 
than the magnetic one for T < 1000 K, and weakly de¬ 
pendent on temperature up to T ^ 500 K. Above this 
temperature, they start to increase, indicating the grad¬ 
ual un-quenching of orbital degrees of freedom. At low 
temperature the frequency dependence of spin- and or¬ 
bital susceptibilities are drastically different. The latter 
shows an activated behavior with a peak at about 0.6 eV, 
that likely comes from transitions from 5" = 1, T = 1 to 
S = 0,L = 2 states, corresponding to an energy differ¬ 
ence 2J. At high temperature, this distinction between 
spin and orbital susceptibilities disappears. 

We also estimated the temperature dependence of the 
entropy. To this aim, we calculated the LDA-I-DMFT to- 
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Figure 2. (a) Active degrees of freedom as revealed by the 
local spin and orbital susceptibilities XmXo (for better vis¬ 
ibility, 10 X Xo is displayed). Inset: inverse susceptibili¬ 
ties (symbols) - there, Xo is multiplied by = 4 for eas¬ 
ier comparison with the high-T local-moment behavior. At 
intermediate T, the spin susceptibility can be fit well by 
Xs = m/{T + Tm) (full line) with Tm = 450 K and m = 2.5fiB 
(close to the atomic S — 1 moment m = S{S+l)g^ (b) 
Frequency-dependence of the local spin and orbital suscepti¬ 
bilities. (c) Total energy (circles; the dashed line is a linear fit 
between 200 K and 800 K) and chemical potential (squares; 
the dashed-dotted line is a 6th order polynomial fit) vs. T. 
(d) Estimated entropy. The arrows indicate the temperatures 
at which the full atomic spin (In 3) and total entropy (spin 
and orbitals. In 9) is reached. 


tal energy as a function of T, displayed in Fig. 2(c) and 
used the thermodynamic relation TdS/dT = dE/dT (see 
supplemental material [30]) to estimate the entropy. It 
is seen (Fig. 2(d)) that the entropy corresponding to un¬ 
quenched spins for Ru^"*" (spin-1, S/ks = In 3) is reached 
by room temperature, while the entropy corresponding 
to both unquenched spins and orbitals (spin-1, L = 1, 
S/kB = In 9) is reached at much higher Torb HOO K. 

These findings that reveal that the decoherence pro¬ 
ceeds in two stages, in which the orbital entropy is re¬ 
leased only after the spins are fully liberated, reinforce 
the Hund’s metal picture of Sr 2 Ru 04 . They are in line 
with the analytical study of Refs [34, 35], in which a 
quantum impurity model appropriate to these systems 
was considered. This model involves three Kondo-like 
coupling constants, corresponding to spin-only, orbital- 
only and mixed spin-orbital degrees of freedom. It was 
found that the spin coupling is much smaller than the 
other two, and can even be ferromagnetic, so that the 
screening of orbital degrees of freedom occurs at higher 
temperature, while that of the spin degrees of freedom is 
controlled by the growth of the spin-orbital coupling and 
occurs at a lower temperature. Incidentally, the quench¬ 
ing of orbital moments at high-T may explain why cal¬ 


culations neglecting the spin-orbit coupling [5, 19] may 
still be accurate for ruthenates down to quite low T, even 
though the bare value of this coupling is ~ 0.1 eV [36]. 

Having identified the relevant degrees of freedom, we 
can now give an interpretation of the thermopower in 
terms of a simple Heikes-like estimate [1, 2]. This was 
previously attempted in Ref. [5] which considered an in¬ 
termediate valence mixing Ru^''"(c?^) and Ru^+(d^), 
for which: an = ^ In + ^Iny^. The authors 
of Ref. [5] observed that the first term gives reasonable 
values for the thermopower provided only spin degrees 
of freedom are retained in evaluating the degeneracies 
^(d"), i.e. T>(d") = 25'„-bl with Sn = 1/2,1,3/2 for d^, 
d'^ and d^, respectively (in agreement with Hund’s rule). 
The problem with this reasoning is that the second term, 
corresponding to conhgurational entropy, is neglected 
and would otherwise yield a diverging result on approach¬ 
ing the actual valence d^ {x = 0) of Sr 2 Ru 04 . We thus 
reconsider the Heikes analysis for this integer valence and 
obtain, by considering the high-T limit of the chemical 
potential (see supplemental material for derivation [30]): 
an = /r/eT = |f In = |f In 2 ~ 30^V/K. This ex¬ 
pression differs from the previous one on several counts. 
Importantly, the problematic configurational term does 
not appear. Furthermore it involves only the degenera¬ 
cies of the two neighbouring configurations d^ and d® 
(note also the factor 1 /2 in front of the logarithm). The 
estimated value ^ 30^V/K agrees reasonably with the 
observed one in the room-temperature plateau regime, 
and points to the heart of the Hund’s coupling dominated 
nature of ruthenates: fluctuating spins and quenched or¬ 
bital moments. Keeping orbital degrees of freedom and 
full degeneracies D{d^) = 4 (spin-3/2,T = 0), D{d^) = 6 
(spin 1/2, L = 1) would instead yield a negative value 
fc_B/2eln4/6 = —17.5pV/K. Indeed a{T) decreases at 
higher T, although this negative value is not reached in 
our calculation perhaps because other atomic multiplets 
become relevant, which are not included in the t 2 g de¬ 
scription. 

Although the agreement between the in-plane ther¬ 
mopower and the entropic Kelvin formula is remarkable, 
one has to keep in mind that this formula is only ap¬ 
proximate and derived in a way that largely ignores the 
transport nature of thermoelectric coefficients. For in¬ 
stance, in a non-cubic system, these coefficients can be 
different along different crystal axis, as recently discussed 
for organic conductors in Ref. [37]. It was pointed out 
there (see also [8]) that the Kelvin formula, which in¬ 
cludes no information about current matrix elements, 
obviously cannot account for such an anisotropy. We 
have calculated the c-axis thermopower of Sr 2 Ru 04 , and 
found a striking manifestation of this effect. As reported 
on Fig. 3, the c-axis thermopower behaves differently 
from both the in-plane one and from the Kelvin approx¬ 
imation above room-temperature. It becomes remark- 
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Figure 3. Temperature-dependence of the calculated c-axis 
thermopower (triangles), compared to that of the in-plane 
one (squares). The errorbars (not shown) are for the c-axis 
response smaller than 5^V/K. The Kelvin and Heikes esti¬ 
mates are also shown, as a plain line and dotted horizontal 
line, respectively. 

ably large for a metallic system, reaching ~ 100/rV/K at 
T ~ 1000 K. Since, to our knowledge, the c-axis ther¬ 
mopower has not yet been measured for Sr 2 Ru 04 , this 
finding is a prediction for future experiments. 

To get basic insight into when the Kelvin formula 
can be expected to work, it is instructive [8] to con¬ 
sider a non-interacting system in which d^/dT\n = 
—ks J deD{e)(e — n)df /de/ f deD{e)df /de with £)(e) the 
density-of-states. The expression of the thermopower has 
the same form, with the important difference that the 
transport function along the appropriate axis replaces 
D{e). Hence, the Kelvin formula is expected to be a 
good approximation only when the transport function 
has a similar energy-dependence as the density of states. 

In Fig. 4 we display the LDA-I-DMFT (a) and the 
LDA(b) transport functions. In contrast to the in-plane 
one, the c-axis transport function exhibits a pronounced 
peak at negative energies ('^ —0.8 eV in LDA, renormal¬ 
ized down to ^ —0.35eV in LDA-I-DMFT). The LDA 
density-of-states (dashed) has a three-peak structure, 
with the outer two peaks originating from the quasi Id 
xz/yz bands. The c-axis transport function exhibits a 
pronounced negative frequency peak only and is thus 
characterized by a strong asymmetry favouring holes. 
This occurs because, as illustrated on panels (b,c) of 
Fig. 4, the dispersion of the bands as a function of kz 
is much stronger close to the center of the Brillouin zone, 
which corresponds to occupied states below Fermi level. 
The large c-axis thermopower is thus due to an inter¬ 
layer ‘hole-filtering’ mechanism. The origin of such a be¬ 
havior can be traced to the body centered crystal struc¬ 
ture, in which the main hopping along the c-direction 
proceeds via the central atom. In a tight-binding pic¬ 
ture, this gives rise to a band energy of the form ej: ~ 
e{kx,ky) — 2tz cos(fcfc/2) cos(fey/2) cos(fcz/2) and the cor¬ 
responding velocities in the z—direction are indeed max¬ 
imal at the zone center where the much larger in-plane 
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Figure 4. (a) Transport functions calculated with 

LDA+DMFT for different temperatures, (b) LDA transport 
function, (c) The band dispersion for kz = 0 and kz = 'k jc as 
a function of in-plane wave-vector. Note that the difference 
between the two energies A is maximal close to kx = ky = 0, 
which corresponds to energies in the range —0.8 to —1 eV. 
This explains the strong peak observed in (c) in this energy 
range, (d) A color intensity map of A in the in-plane Bril¬ 
louin zone, for the band denoted by an arrow in (b). The 
Fermi surface sheets are indicated with lines. 


term e{kx, ky) (counted from Fermi level) is strongly neg¬ 
ative. The effect is largest when the thermal width of the 
function (e — fj,)df/de becomes comparable to the en¬ 
ergy of the peak in the LDA-I-DMFT spectral function: 
bksT ~ 0.35 eV, i.e at temperatures > 800 K. 

In summary, we have shown that the thermopower of 
ruthenates carries key physical information, and indi¬ 
cates in particular that the decoherence temperature at 
which the corresponding entropy is unquenched is much 
smaller for spin than for orbital degrees of freedom, a 
characteristic signature of a Hund’s metal. Our calcula¬ 
tions also predict a different behaviour and a large value 
of the c-axis thermopower at high temperature, which 
results from an interlayer ‘hole-filtering’ inherent to the 
crystal structure of Sr 2 Ru 04 . In contrast to the in-plane 
results, this illustrates the limitation of entropic interpre¬ 
tations of the thermopower. That entropy limitations can 
be overcome may be good news for applications, and it 
would be worth exploring whether such hole-filtering can 
be found and exploited in thermoelectric materials. 
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GENERALIZED HEIKES EORMULA FOR 
INTEGER OCGUPANCIES IN MULTIORBITAL 
SYSTEMS 


in which cIn is the degeneracy of the atomic state with 
N electrons and Aat its energy. Solving this equation for 
/i reads: 


The Heikes formula [I] approximates the Seebeck co¬ 
efficient a as 


aK,aH = -(^ ( 2 ) 

In this expression, /r is the chemical potential and ()at 
indicates that /x/T is to be evaluated in the atomic limit. 
It is actually simpler to consider /x/Tdirectly (instead 
of evaluating it via the thermodynamic relation /x/T = 
— {dS/dn)E and evaluating the entropy from the number 
of possible configurations, as originally done in Ref. [1]). 

If the average occupancy n is non-integer N < n < N+ 
1 (for N integer), /x/T in the atomic limit is controlled 
by the degeneracies of the atomic state with N and iV -|-1 
electrons, and all other valence states can be neglected. 
The occupancy is given by: 

77 ^ - 

( 3 ) 


pPu = dAf(n - N) a{EM+i-EM) 
dM+i{N -I- I - n) 


( 4 ) 


From this expression, it is clear that jdpL reaches a finite 
value in the high-temperature limit (hence that ^ oc T) 
and that the term involving Ee+i — En provides only a 
subleading correction. As a result, one obtains: 


kn , 

an = — log 
e 


dNin — N) 
dN+i{N -I- I - n) 


( 5 ) 


which is the Heikes formula appropriate for fractional 
occupancy N < n < N + 1, as generalized to account for 
atomic degeneracies in Refs. [2, 3] 

For integer occupancy n = N, the above calculation 
must be modified since both neighboring valence states 
TV ± 1 must be retained, as well as of course N itself. The 
expression oi n = N now reads: 


dN-i{N - i)e-P(EM-i-u{N-i)) + (TV + 


( 6 ) 


Solving for ^/T and neglecting as above the subdominant ANALYTICALLY CONTINUED SELF ENERGIES 
corrections in 1/T now leads to the simple expression: AND COMPARISON TO PADE 


ks . 

otH = log 
2e 


dN-i 

dw+i 


( 7 ) 


We note that this generalization of the Heikes expression 
appropriate to integer valence n = N does not involve the 
configurational entropy term, but only the degeneracies 
of the atomic states with neighboring valence ± I. 


The analytical continuation of imaginary frequency 
data remains challenging even with the most recent fast 
and accurate quantum impurity solvers based in contin¬ 
uous time Quantum Monte-Carlo. For the data shown 
in the main text we used the self energies obtained us¬ 
ing stochastic Maximum Entropy approach. In prac¬ 
tice, we analytically continue auxiliary Green’s function 
G(iuj) = l/(iuj — + fi) for a constant jl and ex- 
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Figure SMI. Analytically continued self energies. 


tract the self energy by inverting the resulting G{uj). To 
estimate the error of our procedure, we set jl to two dif¬ 
ferent values (that should give identical answer for ideal 
analytical continuation procedure on noiseless data): (1) 
the real part of the Matsubara self energy at the lowest 
Matsubara frequency (2) the double counting correction 
{U — 2J){n — 0.5) Ri 5.25, where n is the occupancy of 
the t 2 g shell. The two values of p, are close, the former 
being about 0.2eV larger than the latter. As an indepen¬ 
dent crosscheck, we analytically continued the self ener¬ 
gies also directly using Fade approximants. 

The real frequency self energies for selected tempera¬ 
tures are presented in Fig. SMI. One sees that the two 
choice of p (denoted by Maxent 1 and 2, respectively) 
give very similar results. The data agrees quite well also 
with the self energies obtained with Fade, in particular 
for energy |a;| < 0.5eV. 

On Fig.SM2 we show the calculated Seebeck coeffi¬ 
cients using as an input the different real frequency self 
energies just discussed (stars, crosses and squares for the 
two choices of p and Fade, respectively). Whereas quan¬ 
titatively the differences are certainly not negligible, as 


the calculation of the Seebeck coefficient is very sensitive 
to the low energy particle-hole asymmetry of the imagi¬ 
nary part of the self energy, that is particularly difficult 
to extract accurately, all the data gives similar tempera¬ 
ture dependence, with a maximum at about 450K. 

Overall, the maximum entropy calculations agrees well 
with Fade result. Near 450K, however, the Fade data dis¬ 
plays an abrupt jump. The origin of this jump can be 
traced to abrupt shift of the minimum of |ImE(w) to¬ 
wards the negative frequency side in the Fade data. The 
shift of the minimum that occurs at the a temperature 
where the orbital moments start unquenching (see main 
text), which causes the maximum in the Seebeck coeffi¬ 
cient occurs in the Maximum Entropy data in a milder 
continuous fashion (see Fig. SMS). Fade approximants 
are unable to fit the self energy through the crossover, 
hence discontinous jump. 

For the main text we used an average of the two 
stochastic maxent runs in the calculation of Seebeck co¬ 
efficient and the difference between the two results was 
used as an indicative error bar. 











Figure SM2. Calculation of Seebeck coefficient using the Max¬ 
imum entropy (stars and crosses) and Fade (open squares). 
Other data is described in the main text. 
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Figure SMS. Imaginary part of the self energy for the xz 
orbital for several temperatures. Note the progressive shift of 
the minimum of the |ImE| to the negative frequency. 


CALCULATION OF TOTAL ENERGY AND 
ENTROPY 

The total energy is calculated from the charge self- 
consistent LDA-I-DMFT calculation, as described in 
Ref. [4]. The value of the energy and the errorbars on 
the plot in the main text were obtained from the vari¬ 
ance of the total energy in 20 consecutive iterations of a 
charge self-consistent LDA-I-DMFT loop, starting from a 
converged solution. 

The precision of our data is not sufficient to obtain 
the entropy by direct integration of the thermodynamic 
relation 

TdS/dT = dE/dT ( 8 ) 

hence we estimate the entropy in an aproximate way. At 
low-T, the temperature dependence of the energy follows 
a Fermi-liquid form E = "fT^/2, with 7 = 38mJ/molK^ 
for Sr 2 Ru 04 . (Mass renormalization found in experi¬ 
ment is consistent with the one found in LDA-I-DMFT 
[5]). At higher T, the energy is found to depend approx¬ 
imately linearly on temperature. Hence, we used a fit 
E{T) = 7 TV 2 for T < To and E{T) = jTg/2 + c{T -Tq) 
for T > Tq, with c = 7 T 0 such that dE/dT is continu¬ 
ous. We find c ~ O.TSfcs and Tq ~ 166 K. The entropy 
is then obtained by integrating the thermodynamic re¬ 
lation Eq. 8 , which yields S = 7 T for T < Tq and 
S = 7 T 0 -I- 7 Tolog(T/To) for T > Tq. This estimate 
of the entropy is displayed on Fig. 2(d) of the main text. 
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